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We have introduced a separable pairing force, which was adjusted to reproduce the pairing prop- 
erties of the Gogny force in nuclear matter. This separable pairing force is able to describe in 
relativistic Hartree-Bogoliubov (RHB) calculations the pairing properties in the ground state of fi- 
nite nuclei on almost the same footing as the original Gogny interaction. In this work we investigate 
excited states using the Relativistic Quasiparticle Random Phase Approximation (RQRPA) with 
the same separable pairing force. For consistency the Goldstone modes and the convergence with 
various cutoff parameters in this version of RQRPA are studied. The first excited 2 + states for the 
chain of Sn-isotopes with Z — 50 and the chain of isotones with N = 82 isotones are calculated in 
RQRPA together with the 3~ states of Sn-isotopes. Comparing with experimental data and with 
the results of original Gogny force we find that this simple separable pairing interaction is very 
successful in depicting the pairing properties of vibrational excitations. 

PACS numbers: 21.30.Fe, 21.60.Jz, 24.30.Cz, 24.30.Gd 



I. INTRODUCTION 

At present Density Functional Theory (DFT) based 
on the mean-field concept has been widely used for all 
kinds of quantum mechanical many-body systems. In 
nuclear physics the relativistic mean field theory based 
on DFT has been of great success in describing the prop- 
erties of many nuclei all over the periodic table [l| . Con- 
ventional DFT with a functional E[p] depending only 
on the single particle density p can be applied in nu- 
clear physics practically only in a few doubly closed 
shell nuclei. The vast majority of nuclei and in par- 
ticular those far away from the (3 stability line the in- 
clusion of particle-particle (pp) correlations is essential 
for a quantitative description of structure phenomena. 
In the framework of DFT pairing correlations are taken 
into account in the form of Bogoliubov theory 0, H[ for 
the ground states and in Quasiparticle Random Phase 
Approximation (QRPA) for the excited states. Although 
monopole pairing or density dependent delta pairing in- 
teractions are widely used because of their simplicity, a 
cutoff parameter has to be introduced by hand. In or- 
der to avoid the complicated problem of a pairing cut- 
off the finite range Gogny force has been applied [J, Q . 
Its parameters have been adjusted very carefully in a 
semi-phenomenological way to characteristic properties 
of the microscopic effective interactions and to experi- 
mental data @,[3| ■ Over the years the relativistic Hartree- 
Bogoliubov theory (RHB) [5| with finite Gogny pairing 
force has turned out to be a very successful way to de- 
scribe pairing correlations in nuclei. The price we have to 
pay for the advantages is much more numerical efforts in- 
volved, especially in calculations of deformed nuclei and 
in applications for excited states. 

As presented in Refs. [1, @, we have introduced a new 
separable form of the pairing force for RHB theory. A 
similar ansatz has been used in the pairing channel of 



non-relativistic Skyrme calculations in Refs. mini. The 

parameters of our separable force are adjusted to repro- 
duce the pairing properties of the Gogny force in nuclear 
matter. It preserves translational invariance and it has 
finite range. Applying well known techniques of Talmi 
and Moshinsky [12], [H, HH, HH this pairing interaction 
can be used in relativistic and in non-relativistic Hartree- 
Bogoliubov or Hartree-Fock-Bogoliubov calculation of fi- 
nite nuclei. It avoids the complicated problem of a cutoff 
at large momenta or energies inherent in other zero range 
pairing forces. In Ref. [|[ it has been shown that with this 
force the pairing properties of ground states can be well 
depicted on almost the same footing as with the original 
Gogny pairing interaction. 

For excited sates it is important to combine the RPA 
and the RHB in a consistent way in order to describe 
the excitations in unstable nuclei near drip-line, es- 
pecially when the pairing correlations play a crucial 
role. Recently Ring et al. [H, [TtJ have used time- 
dependent relativistic mean field theory to derive the 
fully self-consistent relativistic random phase (RRPA) 
and relativistic quasiparticle random phase approxima- 
tion (RQRPA). For the pairing channel the finite range 
Gogny force D1S is used. Excited states are calculated in 
a consistent framework using the same density functional. 
It has been shown in several applications that RQRPA 
provides an excellent tool for the description of the mul- 
tipole response of stable as well as of unstable and weakly 
bound nuclei far from stability. These investigations have 
been devoted to low-lying collective excitations [H, 
l2fjj |. to giant resonances [2l|, [2^, [H, [H, [25[ , to spin-isospin 
resonances [13, HH, and to new exotic modes in stable [27| 
and unstable nuclei [H, |H, HJ HU, |1] . 

Of course, in the case of spherical nuclei, the calcu- 
lation of QRPA-matrix elements of the original G ogny 
force is possible and computer codes are available \u[ . 
Although the evaluation of such matrix elements for the 
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new separable force is much faster, nonetheless the appli- 
cation of this force for QRPA-calculations in spherical nu- 
clei does not bring an essential advantage. This is, how- 
ever, no longer true for QRPA-calculations in deformed 
nuclei [33] ■ Here one has to do with several thousands 
two-quasiparticle configurations and several millions of 
matrix elements, in particular in relativistic applications 
where the Dirac sea has to be treated properly [16j. A 
separable force is also of considerable advantage in all 
cases where the RPA-problem cannot be solved by di- 
agonalization, as for instance for energy dependent self 
energies in the treatment of complex configurations by 
particle- vibrational coupling [34(. In this case one has 
to work at fixed energy and to solve the linear response 
equations at fixed energy [35| ■ It is well known |36| that 
the dimension of the coupled linear response equations 
scales with the number of separable terms and not with 
the number of two-quasiparticle configurations. There- 
fore a separable force brings essential advantages in all 
these cases. 

So far, the separable pairing force has been used only 
in static applications [8( and in this case was very suc- 
cessful. It is not clear from the beginning, whether one 
can also reproduce the dynamic properties of full Gogny 
pairing in such a simple way, because, in fact, as shown 
in Fig. 6 of Ref. Q both forces are not fully identical. 
In particular there is the problem of Goldstone modes 
connected with translational symmetry. It is well known 
that these modes depend in a very delicate way on the 
properties of the residual interaction. Only in the case of 
full self-consistency these modes decouple fully from the 
rest of the spectrum. This is particular important for 
isoscalar dipole excitations, where the large strength of 
the spurious translational mode can contaminate the low- 
lying El-spectrum considerably. Translational invariance 
is one of the essential advantages of the new pairing force 
as compared to older separable pairing forces such as 
monopole, quadrupole, or other multipole pairing forces. 
However, the new force is presented as a sum over separa- 
ble terms and translational invariance is strictly fulfilled 
only for an infinite number of separable terms. As it has 
been shown in Ref. 8], in static applications this series 
converges quickly and one needs only 8 separable terms 
to get convergence. It is not clear, wether this number 
is large enough for a proper treatment of the Goldstone 
modes. 

This paper is devoted to an investigation of all these 
open questions. The new separable pairing interaction is 
implemented in the relativistic QRPA program and de- 
tails for the calculation of the new pp-matrix elements are 
presented. In order to test the numerical implementation 
of the RQRPA equation with the separable pairing inter- 
action we study the Goldstone modes and the consistency 
of the method. In addition we investigate the question 
whether the dynamic properties of pairing correlations in 
vibrational excitations can be reproduced with the new 
pairing force. As it is known the first 2 + excited states in 
semi-magic nuclei are very sensitive to the pairing gap. 



Therefore we investigate the isoscalar quadrupole exci- 
tations in Sn-isotopes and in N = 82 isotones in the 
RHB+RQRPA approach with the new pairing force and 
compare the first 2 + states with those obtained with the 
full Gogny pairing force. Furthermore we calculate 3~ 
excitations in Sn-isotopes and investigate the sensitivity 
of the isoscalar octupole states to the pairing properties. 

The paper is arranged as follows. The theoretical 
formalism of RHB+RQRPA with the separable form of 
the pairing interaction is presented in Sec. II. The con- 
sistency of the method as well as the Goldstone spu- 
rious modes are investigated in Sec. III. The isoscalar 
quadrupole in Sn-isotopes and in N = 82 as well as 
the isoscalar octupole states in Sn-isotopes are calculated 
in the RHB+RQRPA approach, which are discussed in 
Sec. IV. Finally we shall give a brief summary in Sec. V. 



II. THEORETICAL FORMALISM 

We start with the Sq channel gap equation in sym- 
metric nuclear matter at various densities, 
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(k\Vj p °\k') =-Gp(k)p(k') 
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is the separable form of the pairing force introduced in 
Ref. Q with a Gaussian ansatz p(k) — e~ a2fe2 . A The two 
parameters G and a are fitted to the density dependence 
of the gap at the Fermi surface A(fcp) in nuclear matter. 
Comparing with the Gogny D1S force Q, we obtain the 
parameter set of G = 728 MeV-fm 3 and a = 0.644 fm. 

The RQRPA is constructed in the canonical single- 
nucleon basis, where the wave functions of the RHB 
model have BCS form (for details see Ref. [lTj])- I n 
these calculations the same interactions are used in the 
RHB calculation for the nuclear ground state and in the 
RQRPA equations for the excited states, as well in the 
particle-hole (ph) as in the pp-channel. Since the inter- 
action in the p/i-channel is identical to earlier calcula- 
tions [l3| , we discuss here only the derivation of the ma- 
trix elements of the separable interaction of Eq. @ used 
in the pp-channel of the RQRPA equation in finite nu- 
clei. First, we transform the separable force Eq. ([2]) from 
momentum space to coordinate space and obtain 

V(r 1 ,r 2 ,r' 1 ,r' 2 ) = -GS(R-K') P(r)P(r>) i(l-P ff ) , 

(3) 

where R = i(ri + r 2 ) and r = ri — r 2 are the center of 
mass and relative coordinates respectively, and P(r) is 
obtained from the Fourier transform of p(k), 
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The term <5(R — R') in Eq. §5§ insures the translational 
invariance. It also shows that this force is not completely 
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separable in coordinate space. However, in the basis of 
harmonic oscillator functions the matrix elements of this 
force can be represented by a sum of separable terms 
which converges quickly (for details see Ref. @). 

In the pairing channel we need the two-particle wave 
functions coupled to angular momentum J and the pro- 
jector s(l — P a ) restricts us to the quantum numbers 
of total spin S = and total orbital angular momen- 
tum A = J. These wave functions are usually expressed 
in terms of the laboratory coordinates ri and r 2 of the 
two particles, while the separable pairing interaction in 
Eq. §5§ is expressed in the center of mass coordinate R 
and the relative coordinate r of the pair. Therefore we 
transform to the center of mass frame by using the well 
known Talmi-Moshinsky brackets [12, 0, E3] in the no- 
tation of Baranger [l5| 
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h,n 2 l 2 ;Xfi) = J2 M^l^NL^l-^fx), (5) 
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where 



= {NL,nl,\\n x l x ,n 2 l 2 , A) (6) 
are the Talmi-Moshinsky brackets with the selection rule 



2N + L + 2n + I = 2m + h + 2n 2 + h 



(7) 



Here we need these brackets only for the case A = J. 
We therefore can express the two-body function with the 
quantum numbers 5 = and A = J in terms of center of 
mass and relative coordinates by the sum 
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where j = y/Tj + 1 and s = \ and R NL (R, b R ), R n i(r, b r ) 
are radial oscillator wave functions for the center of 
mass and relative coordinates with the oscillator pa- 
rameters b R — b/-\/2 and b r = b\[2. Finally we find 
the pairing matrix elements of the interaction V\ 2V2 i = 
(n 1 lxji,n2l2j2\V\n 1 >li>ji'>n2'h'32>)j in Eq. © as a sum 
over the quantum numbers N, L, N' , U , n, I, n', and 
I' in Eq. ([5]). The integration over the center of mass 
coordinates R and R' leads to N — N', L — L'. Further 
restrictions occur through the fact that the sum contains 
integrals over the relative coordinates of the form 



Rni{r, b r )Y lm {f)P(r)d 3 r. 



(9) 



They vanish unless / = and L = J. The quantum 
numbers n and n' are determined by the selection rule 
([7]) and we are left with a single sum of separable terms 
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where 

s \ h 3i J j 

xM™ i2 / Rn (r,b r )P(r)r 2 dr . 
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For a Gaussian ansatz of P(r) in Eq. ([4]) this integral 
can be evaluated analytically and we find 
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where the parameter a — a/b characterizes the width 
of the function P(r) in terms of the oscillator length b 
and n is given by the selection rule ([7]): n — n\ + n 2 + 
i(/i +l 2 -J)-N. The results of the RHB + RQRPA 
model will depend on the choice of the effective RMF 
Lagrangian in the p/i-channel, as well as on the treatment 
of pairing correlations. In this work the NL3 effective 
interaction [3?} is adopted for the RMF Lagrangian. In 
the pairing channel we use a separable form in Eq. ([3]) 
adjusted to the pairing part of the Gogny force D1S and 
compare the results of RQRPA calculations with the full 
Gogny force D1S in the pairing channel. 



III. VERIFICATION OF CONSISTENCY 

In the following investigations we solve the 
RHB+RQRPA equations with this separable pair- 
ing force. The Dirac spinors are expanded in a 
spherical oscillator basis with Np — 18 major oscilla- 
tor shells [38]. This leads to a very large number of 
two-quasiparticle(2(7p) pairs and to a huge dimension 
of the corresponding RQRPA matrix. In the practical 
numerical calculations a cutoff energy has to be adopted 
and 2gp-pairs with an energy larger than this cutoff are 
neglected. In relativistic RPA we have two types of 
2qp-pairs and therefore two cutoff energies: Eq p is the 
maximum value of the 2gp-energy for positive energy 
states, and Ec a is the maximum absolute value of the 
2gp-energy for states with one fully or partially occupied 
state of positive energy and one empty negative-energy 
state in the Dirac sea. 

To test the numerical implementations of the RHB + 
RQRPA equations with the separable form of the pair- 
ing interaction we first check the Goldstone modes [3^ |. 
As it is known, the Goldstone modes (often called spuri- 
ous excitations) connected with symmetry violations in 
the mean field wave function have zero excitation en- 
ergy and decouple from the physical states for RPA or 
QRPA-calculations based on a self-consistent mean field 
solution and using the interaction derived as the second 
derivative of the energy density functional (35j . The im- 
portance of a consistent treatment of pairing correlations 



4 



in QRPA calculations has been demonstrated in the non- 
relativistic [Z(J El[ and in the relativistic [l?], EH frame- 
work. The zero-energy Goldstone modes also provide a 
rigorous check of the consistency. Two kinds of spuri- 
ous states have been investigated: one corresponds to 
the violation of particle number in the monopole reso- 
nance with the quantum numbers J T = + ; another is 
connected with the violation of translational invariance 
and the spurious center of mass motion in the dipole 
resonance with the quantum numbers J n = 1~ . There 
should be no response to the number operator since it 
is a conserved quantity, i.e., the Nambu-Goldstone mode 
associated with the nucleon number conservation should 
have zero excitation energy. It is observed that the spu- 
rious state of the number operator in the nucleus 22 O 
disappears when the pairing interaction is treated con- 
sistently in the RHB and RQRPA with our separable 
pairing force. 

For sufficiently large values of the two cutoff param- 
eters Ec P and Ec a the response to the corresponding 
generator of the broken symmetry should vanish for all 
non- vanishing energies. The investigation of the conver- 
gence of the RQRPA results as a function of these two 
cutoff parameters provides a very sensitive verification 
for the numerical performance of the code. In Fig. []Ja) 
we show how the response to the neutron number oper- 
ator in 120 Sn varies with the energy for various values 
of the cutoff parameter Eq p in the range from 40 — 280 
MeV. Here the choice Ec a = 1800 MeV includes almost 
the entire negative-energy Dirac bound spectrum, which 
is large enough to yield a convergent result in the usual 
RRPA calculations. For E Cp = 180 MeV the Nambu- 
Goldstone + mode converges to < 0.1 MeV. The choice 
of the cutoff parameter Ec a has also a pronounced influ- 
ence on the calculated isoscalar monopole response. In 
Fig. [2b) we show the peak energies of the isoscalar giant 
monopole resonance (ISGMR) in 120 Sn as a function of 
E Ca - It saturates for E Ca > 1300 MeV. 

In the dipole channel a large configuration space is 
necessary to bring the spurious 1~ state to zero excitation 
energy. In FigQ] (c) we illustrate the convergence of the 
energy of the 1~ spurious state in 120 Sn. The excitation 
energy of the spurious state is plotted as a function of the 
energy cutoff parameter Ec P for a fixed value of Ec a = 
1800 MeV. We find that the excitation energy goes to 
zero slightly slower in the case of the separable pairing 
force than that with the full Gogny force. This might be 
explained by slightly different ranges of the two pairing 
forces. 

As a consequence we use in the following calculations 
for the solution of the self-consistent RHB equations the 
values E Cp = 180 MeV and E Ca = 1800 MeV. This 
leads to a dimension in the order of 3500 2gp-pairs for 
the RQPRA matrix. 

As we see from the Eq. (jlOl) the separable pairing in- 
teraction is not fully separable in the spherical harmonic 
oscillator basis. We have a sum over the quantum num- 
ber N characterizing the major shells of the harmonic 
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FIG. 1: (a) The response of the neutron number operator in 
120 Sn for five values of the 2qp cutoff energy parameter Ec P - 
(b) The excitation energy of the ISGMR in 120 Sn as a function 
of the cutoff energy parameter Ec a ■ The red circles and black 
squares correspond to the results with the separable pairing 
force and the Gogny force, respectively, (c) The position of 
the spurious 1~ state in 120 Sn as a function of the 2qp cutoff 
energy parameter Ec P - The notations are the same as in 
(b). (d) The position of the spurious state in 120 Sn as a 
function of the number No of separable terms in Eq. (|10[l . 



oscillator in the center of mass coordinate. In order to 
study the convergence with the parameter Nq 7 we show 
in Fig.QJd) the 1~ spurious state in 120 Sn as a function of 
A*o. We find that for nuclei around the line of /3-stability 
( 120 Sn), N = 5 is already large enough to bring the spu- 
rious 1~ state to zero excitation energy. In the previous 
investigation in Ref. Q it was found that a somewhat 
larger value of No = 8 is needed to obtain convergence 
for the ground state properties. 
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IV. RESULTS AND DISCUSSION 

First we investigate the lowest 2 + excitations in Sn- 
isotopes for which experimental data are available using 
the relativistic parameter set NL3 in the pft-channel. In 
Fig.[2]we plot the E2 excitation energies and the B(E2) f 
values for the chain of Sn-isotopes as a function of the 
mass number. We compare the results obtained using 
the new separable interaction @ with the calculations 
by Ansari [19( using the original Gogny force D1S in the 
pairing channel. The results obtained with the separable 
pairing interaction are slightly larger than those calcu- 
lated with the Gogny force. However, the discrepancy 
stays very small within a few percent, and the behavior 
of the E2 and the B(E2) values along the chain of iso- 
topes is well reproduced. These small deviations can be 
understood by the fact that the one term separable pair- 
ing interaction, as we discussed in Ref. 0, yields in the 
ground states slightly larger pairing gaps and therefore 
stronger pairing fields than those found with the original 
Gogny force. Therefore, as expected, also the effect of the 
pairing fields in the excited states is slightly increased in 
the case of the separable pairing force. This conclusion 
is consistent with the pairing properties of the ground 
states described in the RHB calculations. 

In the top panel of Fig. fJJa) we show the average pair- 
ing gaps for protons obtained with RHB calculations for 
the chain of the N — S2 isotones from |§ 2 Zr to \^Yh. 
We find that the separable pairing interaction describes 
the average paring gaps of finite nuclei on almost the 
same footing as the Gogny force D1S, although the pair- 
ing gaps calculated with the separable force are always 
slightly larger than those with the Gogny force. In the 
middle and bottom panels of Fig. [3] we display the E2 
excitation energies and the B(E2) f values for the chain 
of N = 82 isotones as functions of the proton number. 
Again very similar results are found for the Gogny force 



FIG. 3: (a) The proton average gap of the N = 82 isotones 
from io 2 Zr to 70 2 Yb. Dependence of E2 (middle) (b) and 
B(E2) (bottom) (c) on the number of protons for N = 82 
isotones. The B(E2) values are given in units of [10 4 e 2 fm 4 ]. 



and for its separable form in the pairing channel. Exper- 
imental data for the pairing gaps, the excitation energies 
of the lowest 2 + states, and the corresponding B(E2) | 
values are also plotted in Fig. [3] for the chain of N = 82 
isotones. The experimental values of the pairing gap in 
even-even nuclei are calculated by the odd-even mass dif- 
ference with the three-point formula [43| . It is shown that 
the agreement between theoretical predictions and exper- 
imental data is reasonable, except for the nucleus 140 Ce. 
Our results for the RHB-calculations show that the nu- 
cleus 140 Ce with the charge number Z = 58 has a closed 
sub-shell for protons: the 7rl<77/ 2 state is fully occupied 
and we observe in the single particle spectrum for pro- 
tons a large energy gap at the Fermi surface. Therefore 
a very small pairing gap and a large excitation energy of 
the lowest E2 state were predicted, which is inconsistent 
with the experimental data. 

We also investigate in the RHB+RQRPA approach ex- 
cited octupole states of spherical Sn-isotopes. In Fig. Q] 
we plot the first and second excited J 71 " = 3 _ -state and 
the peak energy of the giant octupole resonance for the 
Sn-isotopes as functions of the neutron number. Strong 
low-lying 3~ states are found also in the Sn-isotopes, 
which are consistent with experimental observations. We 
can also observe that the results for the 3~ states calcu- 
lated with the separable pairing force are very close to 
those obtained with the Gogny force in Ref. [2(| • The cal- 
culations with the separable pairing force yield slightly 
larger values of the first and second excited states than 
those with the Gogny force as we have seen it in the 
case of low-lying quadrupole excitations, while both of 
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FIG. 4: Energies of the first (Ei) and second (E2) pygmy 
states of E3 and the peak energies of the Giant Octupole 
states for Sn- isotopes as functions of neutron number. The 
star line is the experimental data of first E3 state of Sn- 
isotope. 



them give almost the same peak energies of the giant 
octupole resonances. This is due to the fact that pair- 
ing collections have rather little influence on the giant 
resonances, but a strong effect on the low-lying excita- 
tions in semi- mag ic nuclei. In comparison with the exper- 
imental data [44|, the RHB+RQRPA calculations both 
with Gogny and separable pairing forces can well describe 
the first 3~ states. This again illustrates that both the 
Gogny pairing force and its separable approximation de- 
scribe the pairing properties of excited states on almost 
the same footing. 



V. CONCLUSION 

In summary we have presented first results for 
RHB+RQRPA calculations in finite nuclei based on a 
new separable force in the pairing channel. This separa- 
ble force is translational invariant and has finite range. 
It contains 2 parameters which are adjusted to reproduce 
the pairing gap of the Gogny force in nuclear matter. In 
the RHB calculations for finite nuclei the two-body ma- 
trix elements in the pp-channel are evaluated using well 
known techniques of Talmi and Moshinsky. The sepa- 
rable form of the pairing interaction can describe the 
pairing properties of finite nuclei in the ground states 
on almost the same footing as its corresponding pairing 
Gogny interaction [8|. Similar techniques are used for 



the evaluation of the two-body matrix elements in the 
pp-channel for RQRPA calculations. The numerical im- 
plementations of the RQRPA with this separable pairing 
force are verified by checking for the separation of Gold- 
stone modes connected with symmetry violations in the 
mean field solutions. The numerical convergence of the 
RQRPA calculations as a function of various cutoff pa- 
rameters is shown for the nucleus 120 Sn. We presented 
applications to the lowest 2 + states and the correspond- 
ing reduced transition rates in Sn isotopes and N = 82 
isotones. The isoscalar octupole excitations in Sn iso- 
topes are also investigated. We found excellent agree- 
ment of our results in comparison with those obtained 
by using the full Gogny force in the pp-channel and with 
available experimental data. 

Therefore we can conclude that this simple separa- 
ble pairing interaction can also be applied in future ap- 
plications of the RHB+RQRPA approach in nuclei far 
from stability instead of the complicated Gogny force. 
In particular this will allow us to use realistic, finite 
range pairing forces also in cases where the numerical 
complexity forced us up to now to neglect pairing cor- 
relations completely, as for instance in recent investi- 
gations of magnetic dipole modes based the tilted axis 
cranking approach [45|], or to restrict us to very simple 
monopole or zero range forces in the pairing channel, 
as for instance in relativistic QRPA calculations in de- 
formed nuclei [HI, [4(| . There are also many extensions of 
relativistic density functional theory beyond mean field, 
which were so far only possible with rather simple pair- 
ing forces, such as applications using projection [471 ] onto 
subspaces with goo d s y mmetries, generator coordinate 
methods (GCM) [48, 49], or investigation of complex con- 
figurations in the framework of particle-vibrational cou- 
pling (PVC) [H, [13]. All these methods require a more 
realistic description of pairing correlations in the future. 
Investigations in this direction are in progress. 
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